Anomalous isotopic effect on electron-directed reactivity by a 3-jjm midinfrared pulse 
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We have theoretically studied the effect of nuclear mass on electron localization in dissociating 
H2 + and its isotopes subjected to a few-cycle 3-/xm laser pulse. Compared to the isotopic trend in 
the near-infrared regime, our results reveal an inverse isotopic effect in which the degree of electron- 
directed reactivity is even higher for heavier isotopes. With the semi-classical analysis, we find, 
for the first time, the pronounced electron localization is established by the interferences through 
different channels of one- and, more importantly, higher-order photon coupling. Interestingly, due to 
the enhanced high-order above-threshold dissociation of heavier isotopes, the interference maxima 
gradually become in phase with growing mass and ultimately lead to the anomalous isotopic behavior 
of the electron localization. This indicates that the multi-photon coupling channels will play an 
important role in controlling the dissociation of larger molecules with midinfrared pulses. 
PACS numbers: 33.20.Xx, 33.80.Wz, 32.80.Qk, 42.50.Hz 



The electronic motion inside the molecule is of fun- 
damental importance in determining the formation and 
fracture of chemical bonds. For more than two decades, 
many efforts have been done to study the electronic dy- 
namics in laser-matter interactions [IH1], aiming at the 
control over ultrafast reactions. With the recent devel- 
opment of laser techniques and attosecond science [SHZ] , 
it has become feasible to steer the electron localization 
in dissociating molecules with the carrier-envelope phase 
(CEP) stabilized few-cycle laser pulses [SHIT] or the se- 
quential ultraviolet and near-infrared pulses [T^lTH] . The 
asymmetric electron localization in molecules can be un- 
derstood as the quantum interference of the populations 
that are resonantly transferred among at least two elec- 
tronics states of different parity [13 HE] ■ 

More recently, the control of electron-directed reactiv- 
ity in midinfrared laser pulses has been explored to en- 
hance the electron localization probability via the match 
between the duration of the few-cycle pulse and the disso- 
ciation time of the molecule [TTHT9] . To hold the control 
efficiency for further heavier nuclei, one may have to use 
the pulse with longer wavelength [TS] . While this consti- 
tutes an important step towards the control of chemical 
reactions in larger molecules, we are wondering whether 
the physical mechanism responsible for the electron local- 
ization in midinfrared pulses remains the same as that in 
the near-infrared regime [9|. Meanwhile, for extending 
the control scheme to the larger molecules, the influence 
of nuclear mass on the reactions must be taken into ac- 
count. However, it remains unclear how the nuclear mass 
affects the control efficiency when the midinfrared pulses 
are applied, though it has been shown that the control of 
electron localization is weakened by the growing mass in 
the few-cycle near-infrared field [20] • 
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To understand the influence of mass in reactions, the 
use of isotopes is one of the practice means. This is be- 
cause different isotopes of a given element have a similar 
electronic structure and chemical properties, but differ- 
ent masses. Indeed, the isotopic effect has been an impor- 
tant issue in the research of laser-molecule interactions, 
such as high-harmonic generation |21j . nuclear dynamics 
detection [22] and high-order above-threshold dissocia- 
tion (ATD) [23]. 

On the other hand, thanks to the recent advance in 
the generation of midinfrared source, few-cycle CEP- 
stable pulses at 3 fim have become available in labo- 
ratory [2H [23 j providing a potential for controlling as 
well as obtaining deeper insight into the electron local- 
ization process in heavier molecules. In this Letter, we 
theoretically study the asymmetric dissociation of molec- 
ular hydrogen by using a few-cycle 3-/xm pulse and report 
an anomalous isotopic effect on the electron-directed re- 
activity. The dissociation process is analyzed with the 
semi-classical model [25]. For the 3-/xm field, the un- 
derlying mechanism for the electron localization is found 
to be based on the interferences between the one- and 
higher-order photon coupling channels. Due to the en- 
hanced ATD for heavier isotopes [23], the interference 
maxima gradually become in phase with growing nuclear 
mass, resulting in the anomalous isotopic effect where the 
electron localization degree is larger for heavier isotopes. 

In our simulations, we have used a numerical model 
that has been well established for studies of electron lo- 
calization in the laser-driven molecular dissociation jSJ[S]. 
This model solves the time-dependent Schrodinger equa- 
tion (TDSE) for the nuclear wavepacket in the Born- 
Oppenheimer (BO) representation. In the present model, 
we have considered the molecular ions oriented along the 
polarization axis. Since the time scale of the nuclear rota- 
tion is several hundred femtoseconds, two reasonable as- 
sumptions underlie our calculation [27j : (i) the nuclei do 
not have time to rotate during the time of interaction and 
(ii) the nuclear rotation after the pulse is not significant. 
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The molecular dynamics of ti2 + can be well described in 
terms of the two lowest-lying electronic states, the lsa g 
and 2pa u states. The TDSE can be written as (atomic 
units are used throughout unless otherwise indicated) [5] : 
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where \i is the nuclear reduced mass, and V g and V u are 
the BO potential curves of the lsa g and 2pa u states, 
respectively. The coupling term V gu equals E(t) x D gu , 
with D gu the electronic transition dipole between the two 
states. The electric field E(t) of the pulse is expressed as 
E(t) = s exp[— 21n2(t/-r) 2 ] cos(u;t + (j>), where e , r, w, 
and <j) denote the electric field intensity, pulse duration 
(FWHM), frequency, and CEP of the pulse, respectively. 

The time-dependent nuclear wave functions corre- 
sponding to the 'left' (I) and 'right' (r) atoms within 
the molecule are defined as f ljr = (x g ± Xu)/v2- The 
kinetic energy release (KER) spectrum S[, r (Ek) is ob- 
tained by Fourier analysis of the real-space wave function 
with the bound states projected out [28]. The final prob- 
abilities P;. r of the electron being localized on the left 
or right nucleus are calculated by integrating 5z >r (J5jt) 
over the kinetic energy E k . In respect that the prepara- 
tion mechanism of H2 + generally indicates an incoherent 
Franck-Condon (FC) distribution of vibrational states in 
experiments [33], we have averaged the observables over 
the initial vibrational states weighted by FC factors in 
the present calculation. 

We study the isotopic effect on the asymmetric dissoci- 
ation of H 2 + with a three-cycle (FWHM), 3000-nm laser 
pulse with an intensity of 1 x 10 14 W/cm 2 . The asymme- 
try parameter of the final electron localization probability 
is defined as A = (P t - P r )/{Pt + P r ). Figures (T[a)-|l[c) 
reveal the asymmetries after FC averaging as a function 
of CEP in the 3000-nm pulse for H2 + and its isotopes, 
respectively. For comparison, a similar simulation is per- 
formed for a three-cycle, 800-nm pulse with the same 
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pulse intensity. The results are shown in Figs. [TJd)-[TJf). 
First of all, it is obvious that the asymmetry is enlarged 
in the midinfrared field for each isotope, as predicted in 
[171 118] . Then, more interestingly, we can obtain two 
contrastive observations from Fig. [lj (i) The asymmetry 
for the 800-nm pulse displays a CEP dependence that 
decreases with growing mass. This is coincident with 
the results of the previous studies [20]. But curiously, 
the asymmetry for the 3000-nm pulse shows an inverse 
isotopic behavior in which the degree of the asymmetry 
is even higher for heavier isotopes, (ii) The curves of 
the asymmetry for the 800-nm pulse show a regular sine 
form, but they become irregular for the 3000-nm pulse. 

To gain insight into the interaction dynamics under the 
3-jum field, we present the absolute asymmetry as a func- 
tion of the CEP and the kinetic energy in Fig. [2] with the 
same laser parameters as in Figs.[j|a)-[]|c). The absolute 
asymmetry is given by A(E k , </>) = Si(E k ,<fi) - S r (E k , cj>). 
Generally, the asymmetry shows a strong dependence on 
both the CEP and the KER for all isotopes. In detail, we 
see, on one hand, two major "undulate tilted stripes" cor- 
responding to the left-right asymmetry. As the nuclear 
mass increases, the stripes gradually become more and 
more precipitous and they are almost vertical for T 2 + . 
Additionally, as indicated by the horizontal dashed line, 
the asymmetry is almost stop at about 1 eV for H2 + , 
but the strong asymmetry extends to higher energy with 
growing mass. Consequently, after the integral over the 
KER, the degree of the asymmetry turns out to be larger 
for the heavier isotope. On the other hand, we can see 
that the stripes of the spectra contain several peaks. In 
the case of H2 + in Fig. pta), the peaks at different KER 
respond inconsistently to the CEP, leading to the irreg- 
ular asymmetry curve after the KER integral, as shown 
in Fig. [ija). But, for T 2 + in Fig. [2^c), the alignment of 
the peaks are almost in step. Therefore, the asymmetry 
for T2 + in Fig.[l|c) becomes closer to a sine-like curve. 

While the KER distributions give us a well understand- 
ing of the features of the CEP-dependent asymmetry in 
Fig. [T] the physical mechanism for the anomalous iso- 
topic behavior of the electron localization has remained 
unclear: Why are there asymmetry peaks at different 
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FIG. 1. The CEP dependence of the FC-averaged electron 
localization asymmetry of H2 4 " and its isotopes, D2 + and T2 + , 
for the pulse wavelengths of 3000 (upper row) and 800 nm 
(lower row). 



FIG. 2. The asymmetry parameter A{Ek,4 l ) as a function 
of the CEP and the kinetic energy for (a) H2 + , (b) D2 + , and 
(c) T2 + in the 3000-nm pulse. The other laser parameters are 
the same as in Figs. [lja)-[ljc) 
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kinetic energy and why do those peaks gradually be- 
come in step with growing nuclear mass? In order to 
uncover the underlying dynamical process, we analyze 
the molecular dynamics with the semi-classical model in 
terms of the quasi-static states [26 , which are related 
to the left and right localized wave functions ^S^ r via 
"01.2 = [(cos 9 ± sin^)*/ ± (cos6> =p sin6')* r ]/v / 2, with 
29 = tan" 1 [-2V gu /(V u - V g )}. The corresponding quasi- 
static eigenvalues can be obtained by 
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It has been demonstrated that the electron localization 
is established by the passage of the dissociating molecule 
through an intermediate regime where the lase-molecule 
interaction is neither diabatic nor adiabatic |26j . This 
process can be evaluated by the transition probability 



P ts (t) = exp{- 



ir[V u (R) - V g (R)\ 
4ioD gu (R)E (t) 



}, 
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where w is the laser frequency mentioned above, E (t) 
is the envelope of the laser field, and R tracks the ex- 
pectation value of the internuclear distance of the time- 
dependent nuclear wavepacket. Then we choose the 
molecular ions H2" 1 " of the 5th vibrational state (v = 5) 
and D 2 + of the 7th vibrational state (v — 7) as the 
investigation targets. The initial equilibrium distances 
(2.634 and 2.641 a.u.) and the vibrational state energies 
(—0.5608 and —0.5603 a.u.) of them are close to each 
other so that we can focus on the influence of the nuclear 
mass on the dissociation process. 

In Fig. [3j from up to down the left and right columns 
depict the time evolution of the electric field E, the quasi- 
state eigenvalues V\ 2, the internuclear distance i?, the 
transition probability P tSl and the localization asymme- 
try for the 800- and 3000-nm pulses, respectively. Ac- 
cording to the semi-classical model, the electron localiza- 
tion is established though the region where < Pt s < 1, 
as shown by the shadows in the figure. We find that 
the evolutions of Pt s in the 3000-nm pulse [Fig. [3^h)] are 
quite different from those in the 800-nm pulse [Fig.|3jg)]. 
This indicates a distinct physical mechanism responsible 
for the electron localization in the 3000-nm field. 

In the case of the near-infrared pulse (the left col- 
umn of Fig. [3]), the nuclei separate slowly [Fig. l3je)]. 
At the early time of the interaction, the energy gap be- 
tween Vi,2 is large and the dynamics is adiabatic. As the 
interaction continues, the energy gap becomes smaller. 
For H 2 + (thick solid curve), the dynamics becomes non- 
adiabatic when the internuclear distance sequentially 
passes through the three- and one-photon crossings, as 
indicated by the dashed lines in Fig. [3^c). This can be 
understood as that the electron localization is established 
by interference between the dissociative populations that 
are generated through the three- and one-photon cou- 
pling channels. After that, the electric field is no longer 
strong enough to overcome the energy gap and the dy- 
namics becomes adiabatic again. For D 2 + (thin solid 



curve), however, due to the slower nuclear motion, the 
pulse becomes much weaker when the internuclear dis- 
tance passes the coupling region. Meanwhile, the en- 
ergy gap is still large for D 2 + [Fig^c)], so the transi- 
tion probability is near upon zero, resulting in the small 
asymmetry. Therefore, under the mechanism of the near- 
infrared field, the electron localization is more difficult to 
be achieved for heavier molecules. 

However, the underlying dynamics is different in the 
midinfrared field (the left column of Fig. [3]). As shown in 
Fig. ^d), the evolutions of Vi j2 are almost the same for 
both molecular ions. According to the transition proba- 
bility in Fig. j3^h) , the dynamics becomes non-adiabatic 
when the internuclear distances pass about R = 3.9 a.u. 
[corresponding to the nine-photon crossing of 3000 nm, 
indicated by the dashed line in Fig. fsj^f)]. Then the elec- 
tron localizations for both molecules begin to be estab- 
lished at almost the same time. Due to the fast separa- 
tion of the nuclei, the dynamics has been diabatic when 
the pulse is turned off and, consequently, the electron lo- 
calization is frozen. Similarly to the situation in 800-nm 
field, the effect of nuclear mass in the 3000-nm field is 
to slow the nuclear motion. But the midinfrared pulse 
is strong and long enough to overcome the energy gap. 
Moreover, a distinct difference is that the higher-order 
multi-photon crossing is open in the midinfrared field. 
The population is therefore transferred between the ger- 
ade and ungerade states through those multi-photon cou- 
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FIG. 3. The time evolution of the physical quantities of the 
electron localization process in the 800- and 3000-nm fields 
for H2 + (thick curve) and D2 + (thin curve), (a)-(b) The 
electric field, (c)-(d) The quasi-state eigenvalues, (e)-(f) 
The internuclear distance, (g)-(h) The transition probability. 
(i)-(j) The localization asymmetry. The shadows indicate the 
regions where the electron localization is established. 
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pling channels. Thus the pronounced electron localiza- 
tion can still be achieved for heavier isotopes. 

Based on the semi-classical analysis above, next we will 
introduce the dissociation model illustrated in Fig. |4|a) 
and then discuss the reason for the anomalous isotopic ef- 
fect in the 3000-nm field. The dissociation process can be 
understood as follows. During the interaction, the outgo- 
ing wavepacket sequentially passes through the 9oj to lu) 
crossings shown by the text-arrows in Fig.ga). Then, 
the populations are resonantly transferred between the 
ls<jg and 2pa u states via absorption or emission of pho- 
tons at the crossing. Note that for the few-cycle pulses, 
some of the crossings may be "switched off" and the 
dissociation pathways would be changed if the CEP is 
varied. This is because the instantaneous electric field 
may be zero when the wavepacket reaches the crossing. 
Therefore, the molecule could finally dissociate via net- 
absorption of one or more photons at different initial 
coupling channels. For example, the dashed arrow in 
Fig. Wa) depicts one of the possible pathways that start 
at the 7a; coupling channel: The wavepacket first absorbs 
seven photons at the 7a; crossing, then passes "5a;" with- 
out emission of photons but emits three photon at "3a;" , 
and finally reabsorbs one photon at "la;" , leading to the 
net five-photon absorption at the initial crossing of "7a;" . 

In Fig. [4^a,) , the centers of the KER distribution for 
different dissociation pathways have been shown by the 
horizontal dashed or solid lines. We also show the FC- 
averaged KER spectra normalized to the maximum value 
for the Bi2 + dissociation under the 3000-nm pulse with 
the CEP of 0.57T and 1.07T in Fig. |IJc). Both spectra 
range from to 2.5 eV and involve several peaks that 
originate from the different dissociation pathways. The 
change of the peak values with the CEP demonstrates 
that the contributions from different pathways are indeed 
affected by the CEP as mentioned above. 

Note that those horizontal dashed and solid lines in 
Fig.BFa) suggest the dissociative populations on the lsa g 
and the 2pa u states, respectively. Due to the large laser 
bandwidth of the pulse, the KER distributions for the 
pathways of different parity overlap in energy, leading to 
the interferences and thus the peaks of the asymmetry 
at different KER. By changing the CEP of the pulse, the 
contributions from different pathways will be modulated, 
so that the electron localization exhibits the CEP- as well 
as KER-dependent asymmetry shown in Fig. [2j 

For comparison, in Fig.gb) we present the asymmetry 
parameter An as a function of kinetic energy for three 
isotopes, where is the absolute value of a cut of the 
asymmetry A(Ek,<f>) normalized to the maximum asym- 
metry at c/> = 0.57T, as marked by the vertical dashed lines 
in Fig. [2] One can see that all the coupling channels, from 
"9a;" to " W , contribute to the asymmetry at lower KER 
between 0-0.55 eV. But the asymmetry of higher KER 
is mainly contributed from the relative higher-order cou- 
pling channels. According to the semi-classical analysis, 
the strength of the relative lower-order coupling channels 
will be shrunk for heavier isotopes due to slower motion 
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FIG. 4. (a) Illustration of the dissociation model for H2 4 " 
in the 3000-nm field. The vertical lines with arrows reveal 
the coupling channels that are noted by the text-arrows. The 
thin gray curves indicate the BO potentials dressed by net ab- 
sorbed numbers of photons. The blue dashed arrow shows one 
of the possible dissociation pathways, (b) The absolute value 
of a cut of the asymmetry A(Ek,4>) normalized to the max- 
imum asymmetry at the given CEP of <f> = 0.5-7T, as marked 
by the vertical dashed lines in Fig. [2] (c) The KER spectra 
normalized to the maximum value for the H2 + dissociation 
under the 3000-nm pulse with the CEP of 0.57T and 1.0-7T. 



of the wavepacket which reaches the coupling region when 
the electric field is weaker. Relatively, the contributions 
of the higher-order channels will become more impor- 
tant to the dissociation of heavier isotopes. This process 
can also be understood as the enhanced high-order ATD 
[23] . Therefore, the asymmetry is gradually enhanced in 
higher energy regions with growing mass, as indicated by 
the horizontal dashed lines in Fig. gb) and Fig. § 

Furthermore, the phase shift of the asymmetry peaks 
at different KER depends on the relative phase between 
the overlapped populations on the gerade and ungerade 
states 16J. For the dissociation pathways starting at 
different initial crossings, the corresponding dissociative 
wavepackets have the different initial phases and thus 
respond inconsistently to the CEP. Consequently, the 
phase of the CEP-dependent asymmetry shifts with the 
KER, resulting in the tilted asymmetry stripes shown 
in Fig. [2] But due to the relative enhancement of the 
higher-order coupling channels for heavier isotopes, the 
CEP response of the asymmetry at higher KER gradually 
becomes consistent to that at lower KER. Therefore, the 
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alignment of the asymmetry peaks is more vertical (KER- 
independent) for heavier isotopes, ultimately leading to 
the isotopic behavior where the degree of electron local- 
ization increases with growing nuclear mass. 

In conclusion, we have studied the electron localization 
in H2 + and its isotopes by the few-cycle 3-/zm pulse and 
have found an anomalous isotopic effect on the electron- 
directed reactivity. Compared to the situation in near- 
infrared fields, a distinct underlying dynamics in the mid- 
infrared pulses is the opening of the higher-order (more 
than three) multi-photon coupling channels. Pronounced 
electron localization, even for heavier isotopes, can be 
achieved through interferences between the dissociation 
pathways that start at the high-order crossings. As the 



nuclear mass increases, the contribution of the high-order 
coupling channels is enhanced and the interference max- 
ima gradually become in step, resulting in the larger 
dissociation asymmetry for heavier isotopes. This unex- 
pected isotopic behavior has provided new insights into 
the electronic dynamics inside molecules. In the future, 
the high-order multi-photon coupling channels will play 
an important role in the control over electron-directed 
reactivity of larger molecules with midinfrared pulses. 
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